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^ I Abstract 

^ ; ROSAT spectra of 11 supersoft X-ray sources ( RX J0439.8-6809, RX J0513.9-6951, RX 

; J0527.8-6954, CAL 87, GAL 83, IE 0035.4-7230, RX J0048.4-7332, IE 0056.8-7154, RX J0019.8 
M ■ +2156, RX J0925. 7-4758, AG Draconis) are approximated with theoretical spectra obtained 
^ '. in LTE models for the atmospheres of hot white dwarfs with line blanketing. The confidence 
^ I intervals of parameters derived from these approximations Tcs, log g, Nh, and R'^/cP are 
I determined. 

^ ■ The results are compared with predictions for a model with stable/recurrent thermonuclear 

• burning on the white-dwarf surface. 

(N : 

' 1. Introduction 

m ■ 
o ; 

\^ [ Supersoft X-ray sources (hereafter "supersoft sources") are a class of X-ray object distin- 

^ I guished by ROSAT observations (Triimper et al. 1991, Hasinger 1994, Kahabka and van den 
^ ■ Heuvel, 1997). The main characteristic of these objects is a very soft X-ray spectrum that falls 
off near 0.5 - 1 keV. Blackbody approximations to the X-ray spectra yield temperatures from 
10 to 80 eV. As a rule, the luminosities derived from such approximations are very high (~ 10^^ 
^ ', erg/s), and often exceed the Eddington limit for solar-mass objects. The first observations of 
I supersoft sources were obtained by the Einstein Observatory, but ROSAT was the first satellite 
, ^ ■ that was able to distinguish them as a distinct class and detect significant numbers of these 
^ ■ sources. Gurrently, about 60 bright supersoft sources are known, located in our Galaxy, the 
Magellanic Glouds, the Andromeda galaxy, and NGG 55 (Greiner 2000). 

One widely adopted model that can explain the nature of at least classical double supersoft 
sources is that of van den Heuvel et al. (1992), in which supersoft sources are close binary 
systems containing a white dwarf and a more massive subgiant secondary that overfills its 
Roche lobe. If the mass of the latter star is approximately twice the mass of the white dwarf, 
the accretion onto the white dwarf occurs on the thermal time scale for the secondary at a high 
rate (~ 10^"^ Mq / yr) , which gives rise to stable nuclear burning on its surface, leading to the 
observed soft X-ray emission. The theoretical possibility of such a situation was predicted in 
(Paczynski and Zytkov 1977) and subsequently studied in detail by a number of authors (Iben 
1982, Nomoto 1982, Fujimoto 1982, Iben and Tutukov 1996). 

The spectral energy distribution of such a supersoft source should be described by the theo- 
retical spectrum of a white-dwarf model atmosphere with the appropriate effective temperature 
Tes, gravitational acceleration g, and chemical composition A (see Section 3). Atmospheres 
with Teff ~ 10^ — 10^ K radiate more efficiently at energies 0.1 - 0.5 keV than a blackbody, 
so that applying LTE model atmospheres without taking lines into account yielded bolometric 
luminosities for supersoft sources that were below the Eddington limit for a solar-mass object 
(Heise et al. 1994). 
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The X-ray spectra of supersoft sources have been approximated using a variety of theoret- 
ical models, from blackbody models to non-LTE models taking line absorption into account 
(Hartmann et al. 1999). The parameters of most sources have been estimated using LTE 
models without including the effect of line absorption. 

Our goal in the current study was to rereduce series of ROSAT observations of known 
supersoft sources using a unified method, in order to derive the physical characteristics of 
these objects in a uniform way, approximating the observed fluxes using theoretical spectra for 
blanketed (i.e. including the effect of lines) LTE model atmospheres for hot white dwarfs. It is 
obvious that considering blanketing of non-LTE model atmospheres for hot white dwarfs would 
be more realistic, but such computations are appreciably more complex and time-consuming 
(Hubeny and Lanz 1995). In addition, we note that the temperature structures of blanketed 
non-LTE and LTE model atmospheres are very similar and differ substantially from non- 
LTE models without line absorption (Anderson 1990). Therefore, we consider blanketed LTE- 
model atmospheres of hot white dwarfs to be more realistic than non-LTE models without line 
absorption (Anderson 1990). 

Recent high resolution (~ 0.06 A) calibration observations of the source CAL 83 by the 
XMM- Newton satellite (Paerels et al 2001) showed that its spectrum is rich in absorption 
lines. This indicates the photospheric nature of the spectrum and demonstrates that any 
detailed analysis of the X-ray spectra of these objects obtained by the new generation of space 
observatories must make use of the method of synthetic spectra. 

The full designations of the sources are presented in Table 2. We will use shortened desig- 
nations consisting only of the first several symbols. 

In all, we studied ten supersoft sources. We also present data for AG Dra; although its 
spectrum is very soft, the resulting parameters are very uncertain and are not interpreted 
further. The parameters of the two sources RX J0527 and RX J0513 are the first derived using 
white-dwarf model atmospheres. 

2. Data reduction 

The sample of sources corresponds to the list of Greiner (2000). A number of objects were 
excluded due to the absence of satisfactory observations or the exceptional softness of their 
spectra, which led to very large uncertainties in their parameters (the results for AG Dra 
demonstrate what occurs in this case). In the end, we used 11 sources and 13 observations in 
our study. 

The spectra of the objects were extracted from a circular region with a radius large enough 
to include the majority of the photons. Information about the background was extracted from 
a zone having an area several times larger surrounding the source. To increase the signal- 
to-noise ratios and obtain a high correspondence to the energy resolution of the detector, the 
spectral flux densities were recalculated to new energy channels that were broader than the 
instrumental channels (Table 1). 

We estimated the physical parameters of the sources by approximating the observed spectra 
with theoretical spectra derived from computations of blanketed LTE model atmospheres taking 
into account interstellar absorption. The free parameters in the fitting were the column density 
of interstellar hydrogen Nh, the effective temperature of the model atmosphere T^s, and a 
normalization factor accounting for the geometrical decrease in the flux. R'^/(P {R is the source 
radius and d is the distance to the source). The gravitational acceleration log g was fixed at 
a preliminarily determined optimal value (this parameter influences the spectrum much more 
weakly than the others; see Section 4). 
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Table 2 presents the main data on the sources and the parameters of the observations used. 



3. Method for computing the hot white-dwarf model atmospheres 

The model atmospheres were characterized by their effective temperature Teg, surface grav- 
itational acceleration g, and a scaling factor for the abundance of elements heavier than helium 
A. Solar chemical composition corresponds to A=l, and a heavy- element content a factor of 
ten lower corresponds to A =0.1. We assumed that the atmosphere was stationary and static, 
and consisted of uniform, planeparallel layers. 

The model atmosphere was found numerically, and the distribution of atoms and ions in 
terms of their excitation and degrees of ionization were determined by the Saha and Boltz- 
mann equations assuming local thermodynamic equilibrium. We used a modified version of the 
ATLAS5 code of Kurucz (1970). The computations included the 15 most widespread elements 
(H, He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Fe, Ni). In contrast to the original code 
of Kurucz, we considered all ionization states. Data on the opacity cross sections were taken 
from (Verner et al. 1996). We also took into account Thomson scattering on free electrons 
(Compton effects in the atmospheres ot white dwarfs weaken as the gravitational acceleration 
increases, and do not influence the continuum spectrum (Madej 1994). We included about 
1200 of the strongest spectral multiplets and lines (for the ionization states corre- sponding to 
the conditions in the atmosphere) selected from the list (Verner et al. 1996), whose effect was 
determined via direct integration. We took the line profiles to be Voigt profiles broadened by 
natural damping, the Doppler effect (taking the microturbulent velocity to be 1 km/s), and the 
Stark effect. The Stark half- width was computed using the simple approximation formula of 
Kurucz and Furenlid (1979). We used a grid of ~ 8000 points in frequency (from 3 TO^^ to 10^^ 
Hz) and 100 points in depth to solve the radiative-transfer equation using the method of Auer 
(1976) for three angles. The upper boundary condition assumed the absence of a rising flux 
at the first surface point. We used a diffusion approximation as a lower boundary condition. 
The temperature structure of the atmosphere was computed using a temperature-correction 
method. 

In some high-lying points in the atmosphere (which exert virtually no influence on the 
emergent flux), the radiation pressure can exceed the gravitational force. Physically, this 
implies the presence of a stellar wind from the surface of the white dwarf. However, the study 
of this wind falls outside the framework of the problem at hand, and we will ignore its effect. 

4. Results of the computations 

We computed a grid of model white-dwarf atmospheres with solar chemical composition 
(Anders and Grevesse 1989) and with effective temperatures 5 -10^ — 1.3 • 10^ K in steps of 10"^ 
K and surface gravities log g =7.0 - 9.5 in steps of 0.5 for comparison with the observations. We 
computed only the models that did not exceed the Eddington limit, i.e., with log g > log gs- 
This limiting surface gravity log gE is specified by the relation 

log ^B = 4.88 + 41og TS, (1) 

where T5 = T^s/W^ K. 

The results of the computations are illustrated in Fig. 1. Figure la shows the temperature 
structures of the blanketed and unblanketed models. Figure lb shows the difference between 
the model spectra without lines (dashed) and with lines (solid). This difference is more clearly 



visible in Fig. Ic, which shows the spectra averaged over intervals of 10 eV. The flux of 
the blanketed model is lower at high energies and higher at softer energies than that of the 
unblanketed model, reflecting the differences in their temperature structures. There are also 
local dips associated with the presence of large numbers of strong lines in these regions. Note 
that the uncertainty associated with our poor knowledge of the magnitude of the Stark effect 
for the spectral lines does not significantly influence the resulting spectral energy distribution. 

Figure 2 shows the dependence of the spectra on variations of the model parameters: ef- 
fective temperature (Fig. 2a), surface gravity (Fig. 2b), and heavy- element abundance (Fig. 
2c). At energies above 0.5 keV, variations in the slope and overall flux of the spectra are most 
appreciable in the presence of variations of the effective temperature. Variations of log g in- 
fluence the overall redistribution of energy much more weakly, and variations of the chemical 
composition are manifest only via differing values for absorption jumps. 

5. Fitting technique 

The PSPC instrument on ROSAT has the following characteristics at soft energies. The 
energy resolution of the detector (AE/E) to 0.5 keV is about film shielding the instrument. 
In addition, absorption 90%. The detector has no sensitivity near 0.4 keV in the interstellar 
medium exerts a large influence at since photons with this energy are absorbed in the these 
energies. If there are no photons with energies higher than 0.5 keV in an observed spectrum, 
variations in the interstellar absorption may be compensated by the joint influence of the 
normalization for the dilution of the flux and the effective temperature, which affects the 
slope of the spectrum. As a result, various combinations of these three parameters can yield 
statistically similar results. This is usually visible in a plot of confidence-level contours in the 
Nh - Teff plane (curves bounding zones within which the true parameter value is located with 
some probability) as a characteristic region of possible T^s and Nh values (Fig. 3). As a rule, 
the upper left part of this zone is consistent with the Galactic value of Nh in the specified 
direction (we will call this the Galactic Nh value) and yields plausible luminosities for the 
source. 

In connection with this, we estimated the errors in T^s and the normalization when limiting 
the column density Nh in the region of most plausible parameter values. We adopted the 
extreme values of the errors of Tes and the normalization for two fixed values of A^^ for the 
boundary values delimiting this region. We first determined the 90% probability boundaries 
for Nh- Then, we determined the errors in the two other parameters for two values of Nh - 
the minimum value of Nh and either the maximum value of Nh or double best fit, depending 
on which turned out to be closer to the best approximation. Observations for which errors 
were derived in this way are denoted in Tables 3 and 4 by the letters A and B in the Bype of 
errors column. Statistically, these errors yield a probability of about 98% that the true value 
is within the error interval. 

For several sources, the parameters were well constrained in the zone corresponding to the 
best-fit approximation. This is because the data for these sources were fairly hard, so that either 
the cutoff in their spectra is beyond 0.5 keV (RX J0513) or the maximum observed flux occurs 
at 0.7 — 0.9 keV (CAL 87 and RX J0925). In this case, we did not use the error-estimation 
method described above, and the parameter errors were calculated using the usual statistical 
criteria (x^ =2.71, 90% probability, type C in Tables 3 and 4). 

Given the complexity of accurately estimating parameters, we attempted to determine the 
most general bounds for the possible parameter values. Most often, due to the impossibility 
of separating the influences of individual parameters, the magnitude of remained virtually 
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unchanged for any value of the surface gravity; however, different effective temperatures were 
naturally obtained for different surface gravities. Therefore, we present two approximations 
for nearly all the sources, corresponding to the lower and upper bounding values of the surface 
gravity. 

In addition, very often, very broad limits for the column density Nh were allowed statisti- 
cally. Our analysis (Section 6) showed that the Galactic column density (Dickey and Lockman 
1990) is located within the confidence intervals for most of the sources (located in the Large 
and Small Magellanic Clouds). Accordingly, we have included approximations for the source 
parameters obtained by fixing the column density to be the Galactic value (the errors were 
calculated in the usual way, error type C). 

6. Discussion 

The approximation results are shown in Tables 3 and 4, which contain the approximation 
data for log g =7.5 and 9.5, and Tables 5 and 6, which contain the data for these surface 
gravities and the Galactic value of Nh- Figure 3 depicts the spectra and confidence contours 
for the surface gravities used in the analyses for two sources. 

Figure 4 compares our results with the data of other authors. Figure 4a shows such a 
comparison for the effective temperatures. Our data for the minimum log g are shown for all 
the sources; for RX J0439, RX J0527, CAL 83, and IE 0035, we have used the approximation 
with the Galactic Nh value, while Nh has been left as a free parameter for the other sources. 
Overall, the agreement is good, though we should point out a number of features. 

First, the temperatures for RX J0527 and RX J0513 are approximately 100 000 K higher 
than those determined by other authors based on approximations of their spectra using Planck 
functions. Second, for the limiting values of T^s for the two hardest sources CAL 87 and RX 
J0925, we used temperatures derived from non-LTE models (lower bounds) and unblanketed 
LTE models (upper bounds) (Hartmann et al 1999, Hartmann and Heise 1997). We can see 
that the non-LTE values are lower and the LTE values (for the non blanketed models) higher 
than our values. The remaining sources approximated by non blanketed LTE models also 
display somewhat higher temperatures than those we derived. Consequently, we conclude that 
including the effect of lines in the LTE models for the white-dwarf atmospheres leads to lower 
temperatures for the supersoft sources. 

Figure 4b compares our values for Nh in the directions of the sources with previously 
published values. We adopted the extreme values from the literature data as the limiting 
values of Nh- If they were known, we used independent values of Nh derived from optical and 
ultraviolet observations. The figure shows our resulting values for the minimum log g. 

We used the data of the Greiner's catalog (2000), as well as from the following papers: IE 
0056 - (Kahabka et al. 1994) {Nh), RX J0048 - (Kahabka et al. 1994), IE 0035 - (Kahabka et 
al. 1999), RX J0019 - (Beuermann et al 1996), RX J0925 - (Motch et al. 1994) {Nh). 

The gravitational acceleration is poorly constrained. We were able to fix it well only for 
the three sources with the highest temperatures: RX J0925, CAL 87 (the only possible value 
for them is log g = 9.5), and RX J0513 (log g =8.4+0.04, errors shown are). The errors in 
their log g values in Fig. 5 were estimated to be ± 0.5. For the remaining sources, values of 
log g from the entire range of parameters for our model grid are statistically allowed. However, 
the model of van den Heuvel et al. (1992) predicts a rigorous relation between T^s and log g, 
according to which stable hydrogen burning is possible without an appreciable increase in the 
white-dwarf radius only within a rather narrow range of accretion rates, called the stable- 
burning strip (SBS) (Nomoto 1982, Fujimoto 1982, Iben and Tutukov 1996). The edges of 



the SBS depend on the mass of the white dwarf. When the accretion rate M is lower, matter 
will accumulate on the surface of the white dwarf, with subsequent explosive burning and the 
ejection of the accreting envelope (novae, recurrent novae, and symbiotic novae). At higher 
values of M, there will be continuous stable hydrogen burning with a luminosity close to the 
Eddington luminosity. The excess matter will be blown out by the optically thick wind from 
the white-dwarf surface, leading to an increase in the effective radius of the photosphere and 
a decrease in the effective temperature (Kato 1996). Using the theoretical relation between 
the mass and radius of a white dwarf (Popham and Narayan 1995), we can reflect the SBS 
from the M — M (Fujimoto 1982) to the T^s — log g (van Teeseling et al. 1996) plane. The 
SBS, Eddington limit, and positions of the sources in these coordinates are shown in Fig. 3. 
The three sources with well-defined values of log g lie inside the SBS (within the errors in the 
parameters). The temperatures of the remaining sources are also consistent with the SBS for 
the case of the minimum allowed surface gravity (log 51 = 7.5 - 8). 

The computations of models for white dwarfs with stable surface thermonuclear burning 
predict an increase in the photospheric radius in the SBS by a factor of two to three compared 
with the radius of a cool white dwarf (Iben 1982, Fujimoto 1982). However, using source sizes 
exceeding the radius of a cool white dwarf of the given mass by a factor of two when reflecting 
the stable- burning strip onto the Teg — log g plane only leads to a shift of the strip downward 
to the left along the strip itself and the Eddington limit. Therefore, in this case, the conclusions 
drawn above that the sources are located in the stable-burning strip remain valid. 

7. Conclusion 

We have carried out an analysis of archival ROSAT observations of 11 known supersoft 
X-ray sources. We have derived the atmospheric parameters T^s and log g for these sources 
by approximating their spectra using computed theoretical spectra for blanketed LTE models 
of hot white-dwarf atmospheres. The resulting parameter values are in agreement with values 
published previously. 

In the next our paper (Suleimanov and Ibragimov 2003), we will carry out an analysis 
of the parameters obtained from the point of view of their consistency with the theory of 
stable/recurrent burning; luminosity, mass, and radius estimates for seven sources; and a dis- 
cussion of the effective temperature - mass relation. 
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Table 1: Grouping of counts in broad channels (from one-third to one-half the energy resolution 
of the detector) used in the analysis. 



Channel 


Energy 


Channel 


Energy 


number 


band, keV 


number 


band, keV 


1 


0.16 - 0.20 


8 


0.87-1.04 


2 


0.21 - 0.26 


9 


1.05 - 1.24 


3 


0.27 - 0.34 


10 


1.25 - 1.46 


4 


0.35 - 0.44 


11 


1.47 - 1.69 


5 


0.45 - 0.56 


12 


1.70 - 1.93 


6 


0.57 - 0.70 


13 


1.94 - 2.23 


7 


0.71 - 0.86 







Table 2: Data on sources and observations used 



Object 


Observation 


a 


b 


Count 


Off-axis 










rate, phot/s 


angle 


RX J0439.8-6809 


rp400161n00 


04'*39"*49^6 


-68°09'01" 


1.40 


1.'52 


RX J0513.9-6951 


rp900398a02 


05'*13"'48^8 


-69°52'00" 


1.947 


0.57 


RX J0527.8-6954 


rp400148n00 


05^27^48^5 


-69°54'02" 


0.1171 


0.061 


CAL 87 


rp400012n00 


05^46™45"0 


-71°08'54" 


0.121 


0.1 


CAL 83 


rpllOlSOnOO 


05^43™33"5 


-68°22'23" 


0.53 


40.29 


IE 0035.4-7230 


rp400299n00 


00'*37™19f0 


-72"14'14" 


0.52 


0.167 




rp400149n00 






0.41 


0.167 


RX J0048.4-7332 


rp600196a01 


00^48™20;^0 


-73°31'55" 


0.17 


20.87 


IE 0056.8-7154 


rp600455a02 


00'*58'"37n 


-71°35'48" 


0.325 


18.3 




rp400300a01 






0.373 


0.38 


RX J0019.8+2156 


rp400322n00 


00'*19"50fl 


+21°56'54" 


1.96 


0.103 


RX J0925. 7-4758 


rp900377n00 


09'*25"'42.^0 


-47''58'"00^ 


0.675 


40.99 


AG Draconis 


rp200689n00 


16'^01"'41.n 


+66°48'10" 


1.09 


1.88 
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Table 3: Approximation parameters for the minimum log g. 



Ubject 


AT 




^05 g 


log [K/a) 


b lux 




iype 01 




in20 „^-2 
iU cm 


XIU J\ 


_2 

cm s 




[\j.2-Z keV j 

— 2 —1 

erg cm s 




errors 


RX J0439 


J —--LI n 1 
^■'^-4.67 


2.80^^:^2 


7.5 


_ _ 9Q 

-26.63t||9 


4.32 • IQ-ii 


9.30/8 


A 


RX J0513 


_ „ . _l_n A7 




^-4-0.15 


rtrt ——-1-0 1 n 
zo.0/_o.o8 


1.09 • 10~^° 


19.4/8 


C 


RX J0527 


^'•"-24.25 


„^_i_9 zlQ 


8.0 


-23.431^06° 


1.26 • 10"''' 


1.47/8 


A 


CAL 87 


69.2ii:J 


8-2ol[]:i 


9.0 


-28.20l[!|^ 


1.17 • lO"'^ 


13.5/8 


C 


CAL 83 


11-0-6.46 


,1 ,-^+0.94 

4-64+0.64 


8.0 


-27.59+1.51 


2.62 • 10 


16.1/8 


B 


IE 0035 




. or.+0 49 


8.0 


-28.891°:^° 


1.67 • 10"^^ 


6.69/8 


C 




'^•°'^-1.89 


J ^ ^j_n ^zi 

4.96l°:|^ 


8.0 


-29.291°:^^ 


8.16 • 10"^2 


3.42/8 


C 


RX J0048 


27.6li^| 


O.00_o,48 


7.5 


-24.551^1^ 


3.09 • 10~^ 


8.08/8 


C 


IE 0056 


-1 9 1 +7.00 
-'-^■-'^-9.24 


2.69l°:i 


7.5 


—9^ ni+^-^^ 


1.20 • 10-9 


9.41/8 


A 




o 71+21.29 
' ■'--2.38 


4-OOlJJ^ 


8.0 


-28.8512-02 


5.31 • 10-12 


13.2/8 


B 


RX J0019 


16.7;?3!35 




7.5 


-23.6212-^2 


4.46 • 10"^ 


18.1/8 


A 


RX J0925 


163;^^ 




9.5 


-26.231°:^? 


2.53 • 10-'^ 


8.85/8 


C 


AG Dra 


0.0595lH^i 


o Qq+0.41 

o.u»_2 59 


7.5 


OO QO"l~14.82 
-^0-0^-0.18 


7.67 • 10-13 


4.21/8 


B 



Note: Fluxes have been corrected for interstellar absorption. 

Type of errors: A denotes error limits based on the limiting values of Nn (98%); 

B denotes error limits based on the left boundary of Nh and 2Nh for the best approximation (98%); 

C denotes errors derived from usual criteria (90%). The error for log g for RXJ0513 is Icr; 

the remaining parameters were derived with this value fixed. The data for AG Dra are approximate. 



Table 4: Approximation parameters for the maximum log g. 



Object 






log g 


log {R/df 


Flux 


XV d.o.f. 


Type of 




1020 cm-2 


xlO^ K 


cm s^2 




(0.2-2 keV) 
erg cm~2 s~^ 




errors 


RX J0439 


t- nn+0.06 
0-UU_5 00 


"^•-'-"J-0.88 


9.5 


-26.79I3I1 


5.03 • 10-11 


9.31/8 


A 


RX J0513 




Not presented: one posible value of log g 






RX J0527 


91 Q+8.00 
^-■^■^-18.55 


3.4313:1 


9.5 


-25.10l|i 


6.72 • 10"^ 


1.47/8 


A 


CAL 87 




Not presented: one possible value of log 5 






CAL 83 




0-00_i 08 


9.5 


-28.1ll?i? 


2.16 • IQ-io 


16.3/8 


B 


IE 0035 


4.62t|^f 


r Qq+0.75 

O-y^-0.50 


9.5 




1.44 • IQ-ii 


6.61/8 


C 




q r-,+2.97 




9.5 


— 9Q 70+0-72 


7.23 • 10-12 


3.44/8 


C 


RX J0048 


99 7+5.1 

zz. 2 


^•-'-'-'-0.63 


9.5 


-25.88lJ;g} 


5.21 • 10"^ 


9.35/8 


C 


IE 0056 


-'■"-'•■'--8.07 


3.onlJ:i 


9.5 


9^ 81+2-46 


3.35 • 10-10 


9.82/8 


A 




q 90 + 17.42 

0.zo_2 21 


4.5911II 


9.5 


-29.29l}|^ 


4.26 • 10-12 


13.1/8 


B 


RX J0019 


n 7+6.20 

-^-■-■'-8.76 


q q-|+2.04 

"^■"^J^-0.60 


9.5 


-25.40+2f 


2.39 • 10-9 


18.2/8 


A 


RX J0925 




Not presented: one possible value of log 5 






AG Dra 


9 r-i+3,99 
-■''-'--2.r)l 




9.5 


-2o.9Gi^:;][; 


3.49 ■ IQ-i- 


4.20/8 


A 



Note: See comments to Table 3. 



in 



Table 5: Approximation parameters for the minimum log g and Galactic Nh- 



Object 








log g 


log {R/df 


Flux 




XV d.o.f. 




1020 cm-2 


xlO^ K 




cm 




(0.2-2 keV) 
















erg cm ^ 


s-1 




RX J0439 


5.60 


r, 79+0.12 

^•'^-0.15 




7.5 


ZD.Z^_Q 20 


7.99 • 


IQ- 


-11 


9.30/9 


RX J0513 


7.24 


5.T4l°:°i 




8.4 




1.71 • 


10^ 


-10 


33.8/9 


RX J0527 


6.31 


4.65l?:?i 




8.0 


-29.181°;^ 


6.90 • 


10" 


-12 


2.17/9 


CAL 87 


7.58 




Unsatisfactory approximation 








CAL 83 


6.33 


5.04lgi? 




8.0 


-28.57lHi 


4.59- 


10- 


-11 


17.3/9 


IE 0035 


6.94 


4.641H? 




8.0 


-28.45lH^ 


3.62- 


IQ- 


-11 


7.58/9 




6.94 


^•"^-0.13 




8.0 


-28.56lH^ 


2.83- 


IQ- 


-11 


6.10/9 


RX J0048 


4.24 




Unsatisfactory approximation 








IE 0056 


6.16 


^•^^-0.11 




7.5 


-27.02;0;26 


2.56 • 


10" 


-11 


10.3/9 




6.16 


Q on+0.99? 
"3-^^-0.09 




8.0 


07 Rq+0.36 
Z(.DO_Q gg 


1.85 • 


10" 


-11 


13.3/9 


RX J0019 


4.20 


q c;7+0.24 
^•^'-0.20 




7.5 


97 7-1 +0.20 
^' ■' -'--0.22 


3.61 • 


10" 


-11 


19.5/9 


RX J0925 


55.6 




Unsatisfactory approximation 








AG Dra 


3-08 


^■•-•''-O.Ki 




7.5 


or n+0.16 
-2G.41_o 1^ 


7.79 


10" 


-12 


4.73/9 



Note: See comments to Table 3. 



n 



Table 6: Approximation parameters for the maximum log g and Galactic Nh- 



Object 






log g 


log {R/df 


Flux 


XV d.o.f. 




lO^o cm-2 


xlO^ K 


cm s~^ 




(0.2-2 keV) 
erg cm~^ s^^ 




RX J0439 


5.60 


o no+o.ig 

"^•'-'^-0.20 


9.5 


-26.47l0i^ 


7.77 • 10^" 


9.31/9 


RX J0513 


7.24 


Not presented: one possible value of log 5 




RX J0527 


6.31 


5.6olS:i 


9.5 


-29.591°;?^ 


6.84 • 10-^2 


2.01/9 


CAL 87 


7.58 


Not presented: one possible value of log 5 




CAL 83 


6.33 


6.20lg:i 


9.5 


-29 m+°°9 


4.59 • IQ-^i 


17.1/9 


IE 0035 


6.94 


r «1+0.14 
^•"-'^-0.27 


9.5 


-30.88lo;o6 


3.59 • IQ-^i 


7.95/9 




6.94 


'^•"'-'-0.30 


9.5 


-28.9810:^3 


2.81 • IQ-^i 


6.83/9 


RX J0048 


4.24 


Unsatisfactory approximation 




IE 0056 


6.16 


q 01 +0.16 
"^•"-•^-0.20 


9.5 


97 o'7+0.35 
^' •'^'-0.29 


2.50 • 10"" 


10.3/9 




6.16 


«r;+0.87 
O.DO_g 26 


9.5 


-27.9lir7^ 


1.81 • 10"^^ 


13.3/9 


RX J0019 


4.20 


4.061°:?? 


9.5 


-28.011°:^^ 


3.47 • 10-^1 


19.3/9 


RX J0925 




Not presented: one possible value of log 5 




AG Dra 


3.08 


1 4 1+0.16 


9-5 


z/.-i.;_2.oi 


5.24 • IQ-"^- 


4.22/9 



Note: See comments to Table 3. 
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Captions to figures 



Figure 1. Difference between blanketed and unblanketed white-dwarf model atmospheres: 
(a) temperature structure (m is the column density), (b) spectra, and (c) spectra averaged over 
10 eV. The model parameters are T^g = 500 000 K, log g =8.5, A=l. The solid curve shows 
the model with lines, while the dotted and dashed curves show the model without lines. 

Figure 2. Atmospheric model spectra as functions of the (a) effective temperature (log g=8.5, 
solar chemical composition, T^g = (1— 9)-10^ K), (b) surface gravity (Teff = 5-10^ K, solar chem- 
ical composition, log g=7.5 - 9.5), and (c) chemical composition (Tpg = 5 • 10''^ K, log g=8.5, 
chemical compositions of 0.25, 0.5, and 1 of the solar value). The spectra are averaged over 
intervals of 10 eV. 

Figure 3. Observed X-ray spectra for two of the studied sources together with the best- 
fitting theoretical spectra and the regions of admissible parameters in the Nh — Tcs plane 
corresponding to the surface gravities chosen for the analysis (see Section 6). Values of the 
residuals are presented together with the spectra. The contours bound the 68, 90, and 99% 
probability zones, and the crosses mark the best-fit values. The source CAL 87 has a relatively 
hard spectrum, and its parameters are localized with certainty; CAL 83 has a soft spectrum, 
and the zone of allowed parameters in the Nh — T^s plane is very extended. 

Figure 4. Comparison of derived parameters with the data of other studies (references 
given in the text): (a) effective temperatures and (b) column densities of interstellar hydrogen. 



Figure 5. Position of the studied sources in the Teff — log g plane. The bold solid curves 
show the strip of stable burning. The dotted curve indicates the Eddington limit. 
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Figure 5: 
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